Electronic superlattices in corrugated graphene 
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We theoretically investigate electron transport through corrugated graphene ribbons and show 
how the ribbon curvature leads to an electronic superlattice with a period set by the corrugation 
wave length. Transport through the ribbon depends sensitively on the superlattice band structure 
which, in turn, strongly depends on the geometry of the deformed sheet. In particular, we find that 
for ribbon widths where the transverse level separation is comparable to the the band edge energy, 
a strong current switching occurs as function of an applied backgate voltage. Thus, artificially 
corrugated graphene sheets or ribbons can be used for the study of Dirac fermions in periodic 
potentials. Furthermore, this provides an additional design paradigm for graphene-based electronics. 



A single layer of graphite, known as graphenej was for 
the first time studied experimentally in 2004 [ij]. This 
started a massive interest in graphene, mainly because it 
is a 2D gapless semiconductor with massless "relativis- 
tic" quasiparticles 0, [1, 0, 0, @, 0j Hi • An unusual integer 
quantum Hall effect [9| and a predicted minimal con- 
ductivity (7 = 4e^/(/i7r) are among the manifestations of 
the linear energy dispersion in the vicinity of the Fermi 
energy. In addition to being a tool for studying funda- 
mental physics, graphene is also of interest for device ap- 
plications (e^. transistors, lenses and NEMS resonators 



Proposed graphene-based devices typically rely on ex- 
ternal electrostatic gates for controlling the electronic 
transport. In this Letter, however, we show that the ef- 
fective potential induced in a graphene ribbon placed on 
a corrugated substrate, can strongly alter the transport 
properties of graphene. This effective potential is deter- 
mined by the local curvature IJ] of the ribbon and intro- 
duces an additional design degree of freedom of interest 
for both fundamental studies of graphene and graphene- 
based electronic devices. For illustration we focus here 
on periodically modulated surfaces — superlattices — 
and show how the corresponding band structure can be 
readily probed by conductance measurements. 

We consider a single graphene ribbon placed on a cor- 
rugated surface and biased by a small applied DC voltage 
as shown in Fig. 1. In the tight binding description the 
ribbon is described by 
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Here denote nearest neighbors (a and h atoms) and 
{{i,j)) next nearest neighbors. When the ribbon is de- 
formed, the matrix elements t and t' change. Thus, if 
tij = VppTT describes hopping between atoms i and j on 
a flat graphene sheet, we find [IF 
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Figure 1: Schematic illustration of the graphene ribbon sys- 
tem investigated in this paper. The DC-biased ribbon of 
width W rests on a corrugated substrate surface of length 
L and height profile h(x,y). 

Here Uy is the vector connecting atoms i and j in the 
undeformed lattice while dy is the corresponding vector 
after deformation. The surface normals are denoted by 
rijQ) . While a general deformation involves both bending 
and stretching, we will here restrict our attention to pure 
bending deformations in one direction, i.e. z = h{x) (see 



Fig. 1). We write the new matrix elements t 
and to second order in &ih we find 
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We note that this expression is similar but not identical to 
the one in Ref. ij]. Inserting the new matrix elements 
in Eq. ([1]) and expanding around the Fermi points K 
and K', results in an effective Hamiltonian with two new 
terms Aes and $eff which are of order aKd^h)^ (cq = 
1.42 A is the lattice constant). For electrons near the K 
point the Hamiltonian is 

i/eff = fiVF^ ■ [-iV + Aeff (x)] + eVa + $ofr(x). (2) 

Here a — xa^ + y<Jy and <Jx,{y) ^tre Pauli matrices, vp — 
l(f m/s is the Fermi velocity, and eVo can be produced 
by the action of a back-gate. 

Unless time reversal symmetry is broken, e.g. by a 
magnetic field, the effective vector potential Aoff(x) at 
the K-point, and its time reversed counterpart at the K'- 
point, only contribute to second order in a^{d'^hY and 
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will be ignored in what follows. The second new term 
corresponds to an effective potential, which for graphene 
bent along the "armchair" direction is 
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In addition, because of the vanishing of the first order 
matrix element of Aofr(x), a third new term correspond- 
ing to a local variation in Fcrmi-vclocity should be con- 
sidered. However, for long wavelengths, kao ^ 1, this 
term can be shown to be much smaller than the effective 
potential $cff- 

The form of the effective potential is simplified if 
we take the shape of the ribbon to be h[x,y) — 
Asai(n'Kx/ L). This approximation captures the qualita- 
tive behavior of a general periodic potential and results 
in the expression 

$eff(x) = Eo{A/aQf{haQ)^{l - coshx) , (3) 

where fcg = 2mT/L and Eq w 0.22 eV. Note that 
$cff(a^) is positive definite (repulsive), and its strength 
varies rapidly with ks. Because of this, large amplitudes 
A/ao 3> 1 are neccesary. Hence, it is important that Eq. 
Q is transformed to a coordinate system s = s{x) that 
follows the graphene ribbon. The relation between s{x) 
and X is 
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where E is the elliptic integral of the second kind and 
A ~ mrA/L. The total length of the graphene sheet is 
then s{L). For simplicity we will from here on write x 
rather than s{x) for the coordinate along the sheet. 

For narrow graphene nano-ribbons the choice of trans- 
verse electronic boundary conditions is of great impor- 
tance. They depend on the configuration of carbon atoms 
along the edge as well as the ribbon width. In this 
Letter, we use boundary conditions for a metallic arm- 
chair edge [l^. The wave vector quantization in the 
transverse direction (y) gives fc„ — m:/W. In the absence 
of fields the wave functions satisfy the Dirac equation 
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with = ±hvp^yk'^ + k^ . For n > the eigenspinors 
are 
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which are two-fold degenerate and for n — they are 
(/>±(fc,y) = [l,±sgn(fc),l,±sgn(fc)]^. 
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Figure 2: (Color online) Band structure for an infinitely long 
and 1 /im wide metallic graphene sheet with armchair edges 
subjected to a periodic potential Vb sin(27ra;/A) (Vb ~ 2.5 
meV) (left panel) and the corresponding number of bands 
crossing each energy (right panel). The linear band calculated 
for kn = (bold line) is not distorted by the periodic poten- 
tial, which is a manifestation of the Klein paradox. A pseu- 
dogap is visible for E ~ 0.02 eV and another at _E ~ 0.04 eV 
roughly corresponding to the free electron energies at A; = ttA 
and k — 2n\. The band structure is in this case symmetric 
for positive and negative energies (not shown). 

and non-degenerate. The factors exp(i(^„(fc)) = (fc -|- 
ik„)/ ^/k^ + k^ correspond to the incidence angle of the 
electrons. 

In the presence of a scalar potential that only depends 
on X there can be no band-mixing although positive and 
negative energy solutions belonging to the same band 
may mix. Thus we look for solutions of the form 



Mx,y) = J dk-f+{k) [0+(fc, y) + (fc,y)] e'''^ 

+ J dk^~{k) [0+(fc,y)-C(fc,y)] 
which, together with Eq. [21 leads to the equation 

[-ihvFidxCTx + iknCTy) + V{x)] 7„(x) = e%{x) 



(4) 



for the two-component spinor 7n(a;) — [7,^(2;), 7~ (a;)] 
with V{x) = eVc + $eff- 

Before analyzing the curvature effects we begin with 
a general discussion of the band structure for periodic 



potentials in the Dirac equation (see also Refs. jlSl. Il9l|). 
For this we consider an infinite graphene strip of width W 
subject to the periodic potential V{x) = Vb sin(27ra;/L). 
In this case Eq. ^ can be be solved numerically and 
the resulting band structure in the reduced zone scheme 
is shown in the left panel of Fig. [51 The right panel 
shows the number of conducting channels as a function 
of energy and one sees that the periodic potential pro- 
duces pseudo gaps whenever this number has a minimum. 
It is interesting to note that the gap sizes increase with 
band index n. This is consistent with our knowledge of 
Klein tunnelling since a potential barrier has no effect on 
massless relativistic particles while for finite mass there 
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Figure 3: (Color online) Conductivity (cr — GL/W) for a 
graphene sheet of length L = 1 fim with a periodic poten- 
tial Vb sin(27ra;/A) (Vb ~ 2.5 meV, A — 100 nm) as a func- 
tion of gate voltage and sheet width W. The pseudo gaps 
at Va = 0.02 eV and Vg = 0.04 eV shown in Fig. [2] corre- 
spond to local minima in the conductivity. The inset shows 
the conductivity along the white dashed line {W — 1 /im) cor- 
responding to the band diagram in Fig. (2] The thick line rep- 
resents an average over nearby gate voltages, which removes 
the Fabry-Perot interferences (see text). The non-averaged 
data is shown as the thin (blue) line. 
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Figure 4: Effective potential along a graphene sheet of a shape 
defined by h{x) = Asin(27ra;/A); ^ = 20 nm, A = 20 nm. 

is an effect which increases with mass (2^ . Here we con- 
sider one-dimensional motion along the graphene ribbon 
in bands corresponding to quantized transverse momenta 
kn- In the effective equation for the longitudinal motion, 
Eq. , these transverse momenta produce a mass term 
that is zero for n — and finite and increasing with n for 
n > 0. 

The band structure in Fig. [2] can be probed by trans- 
port measurements. To demonstrate this we adopt a Lan- 
dauer approach together with a transfer matrix method 
and calculate the transmission probabilities assuming co- 
herent and ballistic transport. The conductance is then 
found from G = (4e^//i) J^n ^" where t„ is the transmis- 
sion probability for channel n and the sum over n runs 
over all channels. The transmission probabilities t„ are 
obtained from the transfer matrices r„ = 11™=! T"™^ 
[2H . Here the interval < a; < s{L) is divided into 
N steps each having a constant potential Vm — V{xm), 
Xm = ms{L)/N. The transfer matrices between slices are 
found from the requirement that the wave functions be 
continuous everywhere (current conservation). The reser- 
voirs on the left and right sides are taken to be infinitely 
wide graphene strips. As will be seen, this gives rise to a 



Fabry-Perot like interference pattern in the conductance, 
due to reflections at the reservoir-ribbon interfaces, as 
the back gate voltage is varied. These fluctuations are 
expected to smear out at finite temperatures and in the 
presence of impurities. As pointed out in Ref. [131 it 
is in general not allowed (as in the case of a 2DEG) to 
introduce a general adiabatic widening of the strip to re- 
move these fluctuations. Such a widening will introduce 
a changing structure of the transverse boundary condi- 
tions and a detailed description of the edge geometry is 
necessary. 

The conductivity of a finite system of length L = 1 /ini 
is shown as a function of gate voltage and strip width 
in Fig. [3l For widths W > L ^ X the conductivity 
does not change appreciably as the the ribbon becomes 
wider whereas for narrow strips strong alteration of the 
conductivity occurs. The inset shows the conductivity 
for parameters corresponding to Fig. [2l The conductivity 
minima agree with the predicted pseudo gaps from the 
infinite structure in Fig. [2l 

Now we consider a finite length graphene ribbon placed 
on a corrugated substrate. We will specifically consider 
narrow strips where the transverse energy level spacing 
Ae„ ~ hvp/W is of the same order as the band edge 
energy, i.e. A ~ VF. For illustration we chose a rib- 
bon placed on a sinusoidally shaped substrate, h{x) = 
Asin(27ra;/A) with A = 20 nm and A = 20 nm. This 
leads to an effective potential (see Fig. |4]) of the order of 
10 meV with an effective wavelength of the order of 80 
nm (consistent with the assumption kao <^ 1). Choos- 
ing a width of = 100 nm thus puts us in the desired 
regime. Figure [5] shows the conductance of a graphene 
ribbon with L — 1 /im as a function of the back-gate volt- 
age (thick blue solid curve) calculated using the transfer 
matrix method described above. The conductance of a 
flat ribbon of equivalent length s{L) is also shown for 
comparison (red dashed curve). Both curves have been 
averaged over nearby points to remove spurious interfer- 
ences arising from the abrupt boundary conditions. Non- 
averaged data for the corrugated sheet is shown as the 
thin (blue) line. 

A comparison of the conductances of flat and corru- 
gated graphene sheets shown in Fig.[5]reveals two distinct 
features. Firstly, there is an asymmetry in the conduc- 
tance of a corrugated sheet with respect to positive and 
negative back-gate voltages because the average effective 
potential $eff is strictly positive. This results in a total 
shift AVg = $cff as well as in changes in the detailed 
structure. Secondly, the effect of corrugation is clearly 
seen to strongly alter the conductance due to conduc- 
tance channels being switched on and off. This can be 
traced to the band structure for an infinite system with 
corresponding parameters shown in the right panel of Fig. 
\5\{N.B. the unaffected n = band has been omitted). 
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Figure 5: (Color online) Conductance of a graphene ribbon as a function of gate voltage (left panel). The conductance of a 
ribbon placed on a corrugated surface (thick blue solid curve) is asymmetric with respect to the sign of the gate voltage in 
contrast to when the surface is flat (red dashed curve). The origin of this asymmetry is that the curvature-induced potential 
is not symmetric. Focusing on positive gate bias we can identify points (A-G) in the conductance curve, corresponding to 
switching on and off conductance in specific channels. The corresponding band edges for an infinite ribbon are shown in the 
right panel. For negative Vg similar features appear in the conductance (not labelled). 
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Figure 6: (Color online) Conductance of a short metallic rib- 
bon with (solid, black) and without (dashed, red) corrugation. 

Finally, we have also considered short graphene ribbons 
where L ^ W ^ X. Figure [6] shows the conductance of 
a ribbon with s{L) « 160 nm and therefore only two po- 
tential maxima (double barrier). Again, comparing with 
a flat system of equal length we flnd that the overall 
features in the conductance change in the same qualita- 
tive ways as described above. In this case, however, one 
should be careful not to confuse structure in the con- 
ductance due to the Fabry-Perot like interferences on the 
one hand with structure due to resonant tunneling on the 
other. In this case both effects are of similar importance. 

In conclusion, we have shown that placing graphene on 
an artificially corrugated surface produces an effective lo- 
cal potential for the graphene electrons. This potential, 
which is related to the local curvature, can be tailored to 
significantly alter the transport properties of graphene. 
Specifically we have considered the effect on the electrical 
conductance of periodic potentials and showed how the 
band structure manifest itself in graphene nano-ribbons. 
Such a relation between transport properties and geo- 
metrical configurations may add to the number of design 
degrees of freedom available for constructing graphene 
based electronic devices. It may also provide an alterna- 
tive transduction mechanism in graphene based NEMS. 

We are grateful for stimulating discussions with Leonid 
Gorelik. Funding was provided by the Swedish Founda- 
tion for Strategic Research (SSF). 



* Electronic address: landr e as.isacsson@chalmers^| 
[1] K. S. Novoselov et al. , Science 306, 666 (2004). 
[2] A. K. Ceim and K. S. Novoselov, Nat. Mater. 6, 183 
(2007). 

[3] M. Katsnelson and K. Novoselov, Solid State Commun. 
143, 3 (2007). 

[4] K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 98, 

076602 (2007). 
[5] K. S. Novoselov et al.. Nature 438, 197 (2005). 
[6] J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and 

C. W. J. Beenakker, Phys. Rev. Lett. 96, 246802 (2006). 
[7] Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 

438, 202 (2005). 
[8] A. L. Vasquez de Parga et al., arXiv:cond-mat/ 

0709.0360, (2007). 
[9] V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 

146801 (2005). 
[10] J. S. Bunch et al.. Science 315, 490 (2007). 
[11] Z. Xu, Q.-S. Zheng, and G. Chen, Appl. Phys. Lett. 90, 

223115 (2007). 

[12] V. V. Cheianov and V. I. Fal'ko, Phys. Rev. B 74, 041403 
(2006). 

[13] V. V. Cheianov, V. Fal'ko, and B. L. Ahshuler, Science 

315, 1252 (2007). 
[14] A. H. Castro-Neto and E.-A. Kim arXiv:cond-mat/ 

0702562,(2007). 
[15] More detailed derivations and results for other boundary 

conditions will be given elsewhere. 
[16] W. A. Harrison, Elementary Electronic Structure (World 

Scientific, Singapore, 2004), revised ed. 
[17] L. Brey and H. A. Fertig, Phys. Rev. B 73, 235411 (2006). 
[18] C. L. Roy, B. Mendez, and F. Dominguez-Adame, J. 

Phys. A 27, 3539 (1994). 
[19] B. F. Samsonov, A. A. Pecheritsin, E. O. Pozdeeva, and 

M. L. Glasser, Eur. J. Phys. 24, 435 (2003). 
[20] A. Calogeracos, N. Dombey, Contemp. Phys. 40, 313 

(1999). 

[21] S. Datta, Electronic transport in mesoscopic systems 

(Cambridge University Press, Cambridge, 1995). 
[22] M. Katsnelson, Eur. Phys. J. B 57, 225 (2007). 



